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Abstract 

Background: The genus Dunaliella (Class - Chlorophyceae) is widely studied for its tolerance to extreme habitat 
conditions, physiological aspects and many biotechnological applications, such as a source of carotenoids and 
many other bioactive compounds. Biochemical and molecular characterization is very much essential to fully 
explore the properties and possibilities of the new isolates of Dunaliella. In India, hyper saline lakes and salt pans 
were reported to bloom with Dunaliella spp. However, except for the economically important D. salina, other 
species are rarely characterized taxonomically from India. Present study was conducted to describe Dunaliella strains 
from Indian salinas using a combined morphological, physiological and molecular approach with an aim to have a 
better understanding on the taxonomy and diversity of this genus from India. 

Results: Comparative phenotypic and genetic studies revealed high level of diversity within the Indian Dunaliella 
isolates. Species level identification using morphological characteristics clearly delineated two strains of D. salina 
with considerable /3-carotene content (>20 pg/cell). The variation in 18S rRNA gene size, amplified with MA1-MA2 
primers, ranged between -1800 and -2650 base pairs, and together with the phylogeny based on ITS gene 
sequence provided a pattern, forming five different groups within Indian Dunaliella isolates. Superficial congruency 
was observed between ITS and rbcL gene phylogenetic trees with consistent formation of major clades separating 
Indian isolates into two distinct clusters, one with D. salina and allied strains, and another one with D. viridis and 
allied strains. Further in both the trees, few isolates showed high level of genetic divergence than reported 
previously for Dunaliella spp. This indicates the scope of more numbers of clearly defined/unidentified species/sub- 
species within Indian Dunaliella isolates. 

Conclusion: Present work illustrates Indian Dunaliella strains phenotypically and genetically, and confirms the 
presence of not less than five different species (or sub-species) in Indian saline waters, including D. salina and D. 
viridis. The study emphasizes the need for a combined morphological, physiological and molecular approach in the 
taxonomic studies of Dunaliella. 
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Background 

Dunaliella, the unicellular microalga, is one of the best 
studied organisms in both general and applied phycology 
for its higher tolerance to extreme conditions of salinity, 
light, temperature and pH, as well as for its richness in 
natural carotenoids, glycerol, lipids and many other bio- 
active compounds [1-4]. Dunaliella salina is reported as 
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the most halotolerant photosynthetic eukaryote with a 
remarkable degree of tolerance from 0.5 to 5 M salt con- 
centrations (30-300 ppt) [2]. This genus naturally inhabits 
saline and hypersaline waters and has a cosmopolitan dis- 
tribution [5] and of the 28 species of Dunaliella, 23 are sa- 
line or hypersaline [3,6-9]. 

Many countries, including India, use D. salina for the 
industrial production of ^-carotene with wide range of 
applications [4,10-12]. Apart from D. salina, D. tertio- 
lecta is used in aquaculture, while many other species 
were found promising for the production of biofuel, 
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bioprospecting of antioxidants, bioactive compounds etc. 
[4,6]. Considering the economic importance, most of the 
studies were mainly focused on the taxonomic, physio- 
logical and biotechnological aspects of the halophilic 
species D. salina [5,13-19] (especially from Indian sub- 
continent) and on the marine species D. tertiolecta. But 
similar exclusive or comparative studies are rarely avail- 
able for other species [20-23], probably due to their 
lesser importance and/or limited distribution. 

Typically the taxonomy of Dunaliella anchors on the 
morphological and physiological features of the organ- 
ism. Apart from the general morphology, salinity toler- 
ance and carotenoid (especially /^-carotene) production 
are the two commonly studied physiological attributes of 
Dunaliella, where considerable variations have been 
accounted at inter and intra-species levels [24,25]. Re- 
cently, Borowitzka and Siva [3] have given a detailed ac- 
count of taxonomic revision of the genus Dunaliella 
with special emphasis on saline species bringing more 
clarity in classification. Dunaliella are unique in having 
a thin plasma membrane instead of a rigid cell wall [26] 
and are able to change their cell shape and volume in re- 
sponse to changes in osmolarity and other growth con- 
ditions [17,27-29]. Due to this high plasticity of cell 
morphology, the traditional practice of species differenti- 
ation merely based on light microscopic observations 
becomes difficult and time consuming. Consequently 
many misidentifications arose in the literature which 
brought in controversies and confusions in the taxo- 
nomic organization of the genus Dunaliella [3,5]. 

Molecular taxonomy emerged as a faster and powerful 
tool as it is consistent and independent from environ- 
mental factors and growth stages [30]. It seems to be an 
advanced and reliable device for the characterization 
and differentiation of morphologically plastic organisms. 
Since 1999, molecular characterization has been found 
promising in the taxonomy of Dunaliella [6,29]. Cur- 
rently 18S rRNA gene [5,31], Internal Transcribed Spa- 
cer (ITS) region [14,23,27] and large subunit of the 
ribulose-bisphosphate carboxylase (rbcL) gene [9] are 
being widely used as effective molecular tools in Duna- 
liella characterization and biodiversity studies. Use of 
these molecular markers has resulted in the re- 
designation of many species [8,32]. Nevertheless, the 
confusion regarding the taxonomy still persists due to 
the misidentifications, and will be there until a complete 
revision is made, using an integrated approach with mo- 
lecular phylogeny, supported by morphological and 
physiological attributes. Many researchers opined to 
have such an integrated approach rather than a single 
system of taxonomic identification [3,32]. 

In India, Dunaliella are found in salt pans, saline and 
hypersaline ponds, lakes, pools etc. as a major primary 
producer. Many species including D. salina have been 



reported [33] to form blooms in salt pans. However, 
taxonomic characterization of Indian Dunaliella strains 
based on phenotypic and molecular traits is rarely avail- 
able. In this background, we conducted a study on the 
characterization of Dunaliella strains isolated from the 
Indian salinas, using morphological, physiological and 
molecular tools and have made an attempt to depict the 
best possible description on Indian Dunaliella spp. 
Based on the results obtained, taxonomic position and 
diversity aspects of the Genus Dunaliella from Indian 
salinas are discussed. 

Results and discussion 

Morphological & physiological parameters 

Morphologically all ten strains of the green biflagellate 
chlorophytes, isolated from 7 different locations along 
the Indian coast (Figure 1 & Table 1), were identified as 
Dunaliella (Figure 2) following the revision of the genus 
by Borowitzka and Siva [3]. Of the 10 strains, 9 were iso- 
lated from hypersaline water bodies and 1 (SI 35) was 
marine. Though purified by agar plating, the cultures 
were not axenic. All morphological characteristics of dif- 
ferent geographical Indian isolates of Dunaliella are 
summarized in Table 2. 

High levels of morphological plasticity in cell shape and 
size was observed among all the 10 Dunaliella, strains, 
but a general consistency in cell size was noticed within 
the range given (Table 3) [34,35]. Among the 10 strains, 
S135 (Calicut, marine isolate), S089 (CMFRI old strain iso- 
lated from Chennai) and S147 (Kutch) were considerably 
larger while strain SI 33 (Kutch) was the smallest. 

In salinity tolerance study (0.5 - 4.5 M NaCl), suffi- 
cient growth (approximately 5-20 million cells/ml in 28 
days from an initial cell density of 15-60 thousand cells/ml) 
was obtained for each strain in different salinities with 
optimum growth at 1.5 or 2.5 M salt concentrations 
(growth rate was 0.1 ±0.05 div.d" 1 during exponential 
growth period), emphasizing that all the strains (including 
the marine isolate S135) are halophilic in nature. Beta 
carotene was quantified in all the isolates (Table 2) at 
normal' and stressed growth conditions. Under stress 
(3.5M NaCl, irradiance of 100-150 [imol photons m" 2 s" 1 ) 
higher level of the pigment (23.4 & 22.9 pg/cell) was 
recorded in the 2 Indian strains, S089 and SI 35 respect- 
ively while for the Australian reference strain D. salina 
CS265, it was nearly 36 pg/cell. The 3 strains turned or- 
ange/red at high salinity (Figure 2, c, m & p). Lower quan- 
tities of the pigment (<2 pg/cell) were observed in the 
strain S133 from Kutch and the 2 Goa strains, S122 
and S125. For the remaining strains it was around 2-7 
pg/cell, under stress. 

Among the many listed attributes, cell size, colour, 
stigma and /^-carotene accumulation are the major traits 
used to discriminate carotenogenic Dunaliella spp. like 
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Figure 1 Sampling locations of Dunaliella strains used in the present study along the Indian coast. 
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D. salina and D. bardawil/salina. Red D. salina (espe- 
cially at high salinity) was reported to have significantly 
large cell size than other common strains like D. parva, 
D. viridis and D. tertiolecta [3]. Limited carotenogenic 
capacity also discriminated other strains from D. salina 
where the latter can accumulate >20 pg /?-carotene/cell 
[3]. Coesel et al [36] and Olmos et al [5] obtained 10 
pg/cell of /^-carotene under non-stressful growth condi- 
tions for the two hyper producing strains of D. salina, 
CCAP 19/30 and 19/18 respectively. In the present 
study, morphological and physiological observations of 
the 2 strains, S089 and SI 35, revealed that they are In- 
dian strains of D. salina. Discrimination derived from 
basic morphology (taxonomic key), characterized the 
remaining strains as D. viridis except S125 (D. rninuta?), 
S133 (D. viridis/D. bioculata?) and S147 {Dunaliella sp.) 
(Table 2). Detailed morphology and physiology based 
studies illustrated considerable diversity within the In- 
dian strains of Dunaliella but a little confusion prevailed 
due to overlapping features with more than one reported 
species (like the cell size increase and ^-carotene 
content of strains S121 and SI 18 at higher salinity, two 
stigmata of SI 33 at lower salinity and pyrenoid charac- 
teristics of SI 25) (Tables 2 and 3). In the present investi- 
gation, molecular characterization was used as a tool to 
resolve the confusion. 

Molecular characterization based on 18S rRNA gene size 

Amplification of 18S rRNA gene with primers MAI & 
MA2 from different Dunaliella isolates in the present study 



gave products with size ranging from ~ 1820-2640 bp 
(Table 1 & Figure 3). The banding pattern observed 
among the present isolates was found matching with the 
reported gene sizes of 18S rDNA using MA1-MA2 pri- 
mers [5,29]. The dissimilarity in product size observed 
among the different isolates could be explained based on 
the presence/absence or difference in the size of introns 
across different species of Dunaliella [5,29]. With regard 
to the presence of 3 types of Group I introns in 18S 
rRNA gene [37], Olmos et al, [31,38] designed a set of 
conserved primers (MAI, MA2 & MA3) and a set of 
species specific primers (DSs, DPs, DBs). They used the 
conserved primers (MAI & MA2) for preliminary differ- 
entiation of various known species of Dunaliella based 
on the size of the PCR product. Subsequently, morpho- 
logically identical Dunaliella strains (e.g., D. salina and 
D. bardawil) got discriminated by position and number 
of the introns [5,29]. Based on these reports 18S rDNA 
of D. tertiolecta (-1770 bp) lacks an intron, D. salina 
(-2170 bp) has only one intron at 5' terminus, D. viridis 
(-2495 bp) has one longer intron again at S terminus 
and D. parva and D. bardawil have two introns (-2570 bp) 
one each at 5' and 3' terminus. Other than these strains, D. 
peircei having -2088 bp (one intron at 5' terminus) was 
also reported. 

In the present study, based on the 18S rDNA gene 
size, clear grouping of all the 10 Indian Dunaliella 
strains was possible (Figure 3 & Table 1). Out of the 10 
strains only 1 strain, S147 (Kutch) produced the shortest 
band (-1820 bp) showing similarity to that reported for 



Table 1 Geographical origin and gene sequence accession details of Dunaliella strains studied in the present work 



Groups 


Strain 
code 


Isolated 
from 


Geographic 
co-ordinates 


Month of 
collection 


Salinity of the 
sampled water 


18S rDNA 
product size 




Genebank 
accession No. 


















18S rDNA 


ITS region 


rbcl gene 




CS265 


Dunaliella salina; Reference strain from CSIRO collection of living microalgae, Australia 


2210 bp 


JN807321 


JN797804 


JN797820 


1 


MBTD-CMFRI-S135 


Sea water, Calicut, Kerala (WC) 


11°15' N 75°46' E 


May 2009 


33 ppt 


2230 bp 


JF708161 


JN797802 


JN797818 




MBTD-CMFRI-S089 


Kelambakkom salt pan, 
Chennai, TN (EC) 


Culture maintained in CMFRI phytoplankton culture collection, 
isolated from Chennai salt pan. 


2210 bp 


JF708173 


JN797806 


JN797811 


II 


MBTD-CMFRI-S1 18 


Salt pan, Nellore, AP (EC) 


14°16' N 80°07' E 


March 2009 


300 ppt 


2290 bp 


JN807316 


JN797808 


JN797813 




MBTD-CMFRI-S086 


Salt pan, Tuticorin, TN, (EC) 


08°47' N 78°09' E 


February 2009 


300 ppt 


2290 bp 


JF708169 


JN797805 


JN797810 




MBTD-CMFRI-S121 


Pulicat salt lake, AP (EC) 


13°40' N 80°11' E 


March 2009 


150 ppt 


2250 bp 


JN807317 


JN797809 


JN797814 


III 


MBTD-CMFRI-S1 15 


Kelambakkom salt pan, 
Chennai, TN (EC) 


12°45' N 80°12' E 


March 2009 


380 ppt 


2550 bp 


JN807315 


JN797807 


JN797812 




MBTD-CMFRI-S122 


Salt pan, Ribandar, Goa (WC) 


15°30' N 73°5T E 


May 2009 


280 ppt 


2550 bp 


JN807318 


JN797799 


JN797815 




MBTD-CMFRI-S133 


Salt pan, Kutch, Gujarat (WC) 


23°50' N 69°39' E 


July 2009 


320 ppt 


2530 bp 


JF708183 


JN797801 


JN797817 


IV 


MBTD-CMFRI-S125 


Salt pan, Pilar, Goa (WC) 


15°26' N 73°53' E 


May 2009 


260 ppt 


2640 bp 


JN807319 


JN797800 


JN797816 


V 


MBTD-CMFRI-S147 


Salt pan, Kutch, Gujarat (WC) 


23°50' N 69°39' E 


April 2009 


1 80 ppt 


1820 bp 


JN807320 


JN797803 


JN797819 



NB: For convenience strain codes used in text included only third part of full strain code (e.g., S086). AP, Andhra Pradesh; TN, Tamil Nadu; WC-west coast; EC, east coast. 
Indian strains were grouped into subsets based on the 18S rDNA size obtained by PCR amplification with MA1-MA2 primers. 
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Figure 2 DIC microscopic images of different Dunaliella isolates, (a) Dunaliella sp. S086 (Tuticorin salt pan), (b) & (c) D. salina S089 (CMFRI old 
strain), (d) D. viridisl SI 1 5 (Chennai salt pan), (e) & (f) Dunaliella sp. S1 18 (Nellore salt pan), (g) & (h) Dunaliella sp. S121 (Pulicat lake), (i) D. viridisl 
S122 (Goa salt pan), (j) Dunaliella sp. S125 (Goa salt pan), (k) Dunaliella sp. S133 (Kutch salt pan), (I) & (m) D. salina S135 (Calicut marine isolate), 
(n) Dunaliella sp. S147 (Kutch salt pan), (o) & (p) D. salina CS265 (Australian reference strain). In brackets given the origin of isolates, (c) & (m) 
orange red cells of Indian isolates of D. salina (S089 & S135) grown at 4.5 M NaCI concentration, (f) & (h) large yellow green cells of S1 18 and 
S121 at 4.5 M NaCI. (p) Reference strain D. salina CS265 at 2.5 M NaCI turning orange. Scale bar given - 5 urn. 



D. tertiolecta (-1770 bp) probably due to the absence of 
any introns (Group V). The 2 Indian strains S089 
(CMFRI old strain, Chennai) and S135 (Calicut marine 
isolate) and the reference strain CS265 (D. salina) pro- 
duced -2200 bp size band (Table 1) closer to the 
reported D. salina (-2170 bp). This further supported 
phenotypic identification of the above 2 Indian strains as 



D. salina (Group I). Studies using the 18S PCR products 
revealed a clear separation of morphologically similar 
strains (D. viridis?), into 2 groups - (Group II & Group 
III in Table 1 & Figure 3). The 18S rDNA size (-2300 bp) 
of Group II strains (S086 (Tuticorin), S118 (Nellore) & 
S121 (Pulicat)) was showing an indication that these 
strains are more close to D. salina than D. viridis. While 



Table 2 Morphological and physiological characteristics of 10 Indian Dunaliella strains 



Groups 



Strain 
Code 



Cell 
colour 



Cell 
shape 



Flagella 
length 



Stigma Pyrenoid Refractile Mode of 

granules reproduction 
observed 



P Carotene Salinity Identified as 
normal/stress optimum 
(pg/cell) 



MBTD-CMFRI-S1 35 Green to red Ovoid, spherical, 1 .3 or 1 .5 to 

cylindrical cell length 

MBTD-CMFRI-S089 Green to red Ovoid, spherical 1.3 to cell length 



MBTD-CMFRI-S1 18 Green to 
orange 

MBTD-CMFRI-S121 



MBTD-CMFRI-S086 
MBTD-CMFRI-S1 15 
MBTD-CMFRI-S122 
MBTD-CMFRI-S133 



Green 

Green 

Green 

Green 

Yellow 
green 



IV MBTD-CMFRI-S125 Green 



V MBTD-CMFRI-S147 Green 



Ovoid spherical 

Ovoid pyriform 

Ovoid, oval 
or pyriform 

Ovoid, oval 
or fusiform 

Oval, cylindrical, 
Fusiform, Elliptical 



1.5 to cell length 



Not clearly visible 
or diffuse 



Not visible/ 
Diffuse large 



One; large, red, 
median, diffuse 



1.5-2 to cell length One; large, red, 
median, distinct 

1.5- 2 to cell length One; Small, red, 
median, distinct 



1.3 to cell length 
1.3 to cell length 
1.3 to cell length 



One; small, red, 
anterior, distinct 

One; large, red, 
anterior, distinct, 

One, Two at 
lower salinity; 

small, red, 
median, distinct 



Large with 

distinct 
amylosphere 

Large with 

distinct 
amylosphere 

Small with 
amylosphere 

Large with 
amylosphere 

Small with 
amylosphere 

Small with 
amylosphere 

Large, 
Amylosphere 

Small, with 
amylosphere 



Absent Cell division 



Absent Cell division 



Absent 
Absent 
Present 
Absent 
Present 
Absent 



Cylindrical fusiform Equal or 1.3 to 
cell length 



One, large, red, Small, with distinct Absent 
anterior, distinct separate 
starch grains 



Oval, fusiform 



1.5 or 2 to 
cell length 



One, large, red, 
median distinct 



Large with 
amylosphere 



Present 



Sexual, 
cell division 

Sexual, 
cell division 

Sexual, 
cell division 

Palmella, 
aplanospores 

Palmella stage 



Cell division, 

Palmella, 
aplanospores 

Cell division 



Palmella 
(dominant stage), 
Cell division 



8.68/22.94 1.5M NaCI D. salina 

6.53/23.36 1.5M NaCI D. salina 

2.1 1/3.47 2.5M NaCI Dunaliella sp. 

1 .59/2.1 7 1 .5M NaCI Dunaliella sp. 

2.68/3.41 0.5 M NaCI Dunaliella sp. 

1.05/1.99 1.5MNaCI D. viridis? 

0.67/1.78 1.5MNaCI D. viridis? 

0.51/1.26 1.5MNaCI D.viridis/ 
D. bioculata? 

0.70/1.8 1.5MNaCI D. minuta? 

0.89/6.7 1 .5M NaCI Dunaliella sp. 



Grouping of the subsets was formed based on common morphological features including cell size and /3-carotene accumulation at high salinity and light (stress). 



Table 3 Descriptive statistics of cell size variables and F-values (derived from the analysis of variance) of different Dunaliella isolates from Indian coast 



Groups 








II 






III 




IV 


V 




Strain code 


MBTD-CMFRI-S1 35 


MBTD-CMFRI-S089 


MBTD-CMFRI-S086 


MBTD-CMFRI-S118 


MBTD-CMFRI-S121 


MBTD-CMFRI-S1 15 


MBTD-CMFRI-S1 22 


MBTD-CMFRI-S1 33 


MBTD-CMFRI-S1 25 


MBTD-CMFRI-S1 47 


F value 


Length pim 


17.51 ±1.78 


14.12+2.25 


9.15+1.02 


9.51+1.09 


9.37+1.30 


9.02+0.96 


8.46+1.12 


7.91+0.93 


9.89+1.37 


11.17+1.50 


138.33* 




(12.30-21.17) 


(10.01-18.82) 


(6.44-10.68) 


(7.96-12.25) 


(6.45-11.77) 


(6.79-12.12) 


(5.62-10.55) 


(6.54-9.78) 


(8.38-12.99) 


(8.02-13.83) 




Width [xm 


10.30+1.96 


9.57± 1.35 


6.14+0.92 


6.91+0.74 


5.94+0.96 


5.09+0.77 


4.74+0.48 


3.89+0.60 


4.34+0.69 


7.23+1.15 


125.85* 




(8.61-19.79) 


(7.46-12.58) 


(3.52-8.08) 


(5.84-8.76) 


(4.14-7.54) 


(3.02-6.98) 


(3.91-5.76) 


(3.11-5.10) 


(3.27-5.95) 


(5.40-10.07) 





Measurements are presented as, Mean ± SD (min. - max.); ^Significant at the 1% level; SD is standard deviation. 
Grouping of subsets was statistically formed based on the average length/width of the Dunaliella cells. 
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Figure 3 18S rDNA amplification with MAI & MA2 primers in 

1% Agarose gel. Lane 1-1 1 CS265, MBTD-CMFRI-S089, S135, S086, 

S1 18, SI 21 7 S125, S1 15, S122, S133 and S147 respectively & Lane 

M -250 bp ladder (Genie, India), 
v J 



group III strains (S122 (Goa), S115 (Chennai) & S133 
(Kutch)) gave a band size of -2550/2530 bp which could 
be compared to the reported D. viridis (2495 bp) or D. 
parva (-2570 bp) probably with 1 or 2 introns. The 
band size of the Goa strain S125 (-2640 bp, Group IV) 
was however not in accordance with any of the reported 
species of Dunaliella [5]. Partial (-600 bp) sequencing of 
S terminus region of the PCR products could not con- 
firm the presence of any introns, while the generated 
partial sequence information (refer Table 1 for GenBank 
accessions) was found to be highly conserved across spe- 
cies and therefore could not specify the species delinea- 
tion. Further characterization was carried out based on 
molecular phylogeny of a more variable ITS region and 
a conserved rbcL gene for more clarification about spe- 
cies lineages of Indian Dunaliella. 

ITS phylogeny 

The phylogenetic analysis based on ITS region (-700 bp) 
using maximum likelihood confirmed high level of genetic 
diversity within Indian Dunaliella isolates. All Dunaliella 
spp. (including the sequences from NCBI, Table 4) were 
found to be separated into 3 major clusters, with Chlamy- 
domonas reinhardtii forming an out group as expected 
(Figure 4). Majority of the named species of Dunaliella 
from NCBI were found to be grouped under a single clus- 
ter (clade 1) except D. salina CCAP 19/18, D. salina/D. 
viridis CCAP 19/3, D. viridis CONC 002 and Dunaliella 
spp. ABRIINW Ml/2, SPMO 200-2, SPMO 600-1 and 
the reference strain of D. salina CS265. When out 
grouped with D. tertiolecta (Figure 5), the ITS tree 
branching was found well supporting the morphology and 
18S rRNA gene size based grouping (Group I-V) of the 10 
new isolates of Dunaliella. 

The genetic divergence values observed among clade 2 
(Figure 4) isolates ranged up to 9.1% (between S089 & 
CCAP 19/3), which was comparable to that observed be- 
tween different species of the genus Dunaliella [14,17]. 
The 2 Indian D. salina strains (Group I, Figure 5) S089 



Table 4 Sequence accession no.s of Dunaliella and other 
strains from NCBI database included in present study 



Strain 


ITS Accession 
No. 


rbcL Accession 
No. 


D tertiolecta CCAP 19/27 


EF473748 




D tertiolecta ATCC 30929 


EF473742 




D tertiolecta SAG 13.86 


EF473738 




D tertiolecta Dtsi 


EF473730 




D tertiolecta CCMP 1302 




DQ3 13205 


D tertiolecta UTEX 999 




DQ3 13203 


D peircei/D viridis UTEX 2192 




DQ313196 


D primolecta 


DQ 116745 


DQ 173090 


D primolecta 




AB1 27992 


D primolecta UTEX 1000 




DQ313198 


D parva 


DQ 116746 




D parva/D viridis UTEX 1983 




AJ001877 


Dbioculata UTEX 199 


DQ377086 


DQ313195 


D bioculata 




AB1 27991 


D salina Dsge 


EF473732 
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EF473746 
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DQ377106 
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DQ377120 


DQ313218 


Chlamydomonas reinhardtii 


AB5 11842 




Paulschulzia pseudovolvox 




D86837 



(CMFRI old strain) and SI 35 (Calicut marine isolate) got 
clustered with the Australian D. salina strains CS265 
and CCAP 19/18 with divergence values ranging from 
1.9% (between CS 265 & CCAP 19/18) to 5.6% (between 
S089 & S135). Whereas, the strains S086 (Tuticorin), 
S121 (Pulicat) and SI 18 (Nellore) were found closer to 
D. salina/D. viridis CCAP 19/3. The much higher diver- 
gence (>8%) of the 3 Indian strains from the caroteno- 
genic D. salina strains (CCAP19/18 and CS265) was in 
agreement with the grouping of the 3 strains in Group II 
(Figure 5) based on the morphological, physiological and 
18S rDNA size based analyses. 

The remaining 5 Indian Dunaliella strains (SI 15, 
S122, S125, S133 and S147) along with D. viridis CONC 
002 formed a separate cluster (clade 3 of Figure 4). The 
strains showed divergence range from 0% (between SI 15 
& S122) to 7.6% (between S147 and D. viridis CONC 
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002). The 2 D. viridis(?) strains SI 15 (Chennai) & S122 
(Goa) and the D. viridis/D. biocuata(?) strain S133 
(Kutch) were found in close proximity (mean divergence 
of 2.22%) with CONC 002 D. viridis (Group III of Figure 5). 
The other 2 strains S125 and S147 (Group IV & V of 
Figure 5) were found to be well separated from the 
above group with divergence values of 4.98% and 6.42% 
respectively with D. viridis CONC 002. 

The mean pair wise genetic distance values observed 
among the Indian isolates of the 2 major clades (5.35% 
for clade 2 and 5.12% for clade 3) were comparatively 
higher than that observed among the named species of 
Dunaliella (1.14% clade 1). Further, the genetic diver- 
gence values observed among the Indian Dunaliella 
isolates based on ITS sequence variations were consider- 
ably higher than that reported in Chlamydomonas spp. 
(a minimum of 3.5% between 2 species) by Coleman & 
Mai (1997) [39]. Thus, the pattern of genetic divergence, 
along with the phylogenetic divergence pattern, clearly 
indicates the presence of at least 5 or more number of 
species/sub-species among the 10 Indian strains (includ- 
ing D. salina and D. viridis), 

rbcL gene phylogeny 

The pattern of genetic diversity observed among the In- 
dian Dunaliella strains based on rbcL gene sequence 
variations was in accordance with the above observa- 
tions based on 18S rDNA and ITS analysis except for 
the positioning of S147 and SI 35. The phylogenetic tree 
constructed, using maximum likelihood (Figure 6) analysis 
with rbcL gene sequence data, was forming 2 major 
clusters with Paulschulzia pseudovolvox as out group. 
The mean genetic divergence value observed between 
the 2 clades was 5.89% and that observed among dif- 
ferent isolates of Dunaliella ranged from 0.16% to 
7.73%. 

Being a protein coding gene, the pairwise genetic di- 
vergence (Tamura 3 parameter) values observed among 
Dunaliella isolates based on rbcL gene sequences were 
found to be less in comparison with that observed in 
ITS (a non coding region) sequences. The independent 
phylogenetic analyses using ITS (Figure 4) and rbcL 
gene (Figure 6) sequences (Kimura 2 and Tamura 3 
parameters respectively) were found to be taxonomically 
incongruent especially in clade 1. The major topological 
change observed was the change in the positioning of 
the isolate S147. Within rbcL phylogeny, this strain from 
Kutch was found closely allied with clade 1 (Figure 6), 
whereas, with ITS data it was close to D. viridis CONC 
002 and other Indian isolates (S125, S133, SI 15 and 
S122) in clade 3 (Figure 4). Similarly, the marine D. sal- 
ina strain SI 35, got clustered with non-carotenogenic 
strains in sub-clade C (with divergence of 1.15%) instead 
of carotenogenic D. salina (sub-clade B). 



Clustering of all the remaining 8 strains in both ITS 
and rbcL phylogenies was more or less similar. As 
expected S089 (D. salina, CMFRI old strain) clustered 
with carotenogenic D. salina species CS265 and CCAP 
19/18 in clade 1, (sub-clade B, with 100% similarity). 
The positioning of 3 strains S086, SI 18 & S121 (sub- 
clade C, in Figure 6) along with D. salina/D. viridis 
UTEX 200/ CCAP 19/3 and D. peircei/D. viridis UTEX 
2192 (with <1% divergence value) strongly suggests fur- 
ther taxonomic revision. Similarly, in clade 2, the posi- 
tioning of the Goa isolate S125 (with maximum 
divergence 7.33%) and the clustering of strains SI 33, 
S122 & SI 15 with D. viridis CONC 002 and D. parva/D. 
viridis UTEX 1983 (with divergence values 3.15% & 
0.33% respectively) was in concordance with ITS 
phylogeny. 

Grouping of Indian Dunaliella strains 

The 2 larger carotenogenic strains (>20 pg/cell f$- 
carotene content) S135 and S089 forming the GROUP I, 
produced 18S rDNA size -2200 bp and clustered with 
D. salina CCAP 19/18 and CS265 in ITS phylogeny. 
These results confirmed the taxonomical identity of the 
2 strains as D. salina (Section Dunaliella). But closeness 
of SI 35 to the 2 morphologically dissimilar, lower ^-caro- 
tene content strains, S121 (Pulicat) and S086 (Tuticorin) 
in rbcL phylogeny has to be noted, which may be due to 
its marine origin. 

GROUP II included the strains S086 (Tuticorin), SI 18 
(Nellore) and S121 (Pulicat), which clustered with D. sal- 
ina/D.viridis (CCAP 19/3) in ITS phylogeny and had 
-2300 bp band for 18S rDNA. The present study shows 
the closeness of these 3 strains to D. salina by molecular 
analysis (18S rDNA size and ITS & rbcL phylogenies) ra- 
ther than by morphological features. These strains were 
with lesser ^-carotene content (-2-4 pg/cell) and cells 
were always green (only SI 18 turned slightly orange at 
higher salinity), smaller and with a clear stigma, which 
were not corresponding with hyper /^-carotene producer 
strain of D. salina and are more or less characters of D. 
viridis [3]. However, there is a description of a greener 
D. salina (KCTC10654BP) from Korea [17] with low cel- 
lular ^-carotene. But 18S rDNA size details are not avail- 
able for the above Korean strain for comparison. All 
these factors along with the appearance of D. viridis/D. 
peircei UTEX 2192 close to SI 18 in rbcL phylogeny 
(clade 1, sub-clade C) and 18S intron phylogeny of D. 
peircei UTEX 2192 by Hajezi et al. [29], emphasizes a 
need of revisiting the taxonomic identity of all the above 
reported strains along with the 3 Indian strains using 
molecular approaches. 

GROUP III was formed by 3 strains, SI 15 (Chennai), 
S122 (Goa) and S133 (Kutch) allied to D. viridis. This fur- 
ther confirms the possibility of the former 2 strains to be 
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Figure 4 Phylogenetic tree of the Maximum Likelihood (ML) analysis inferred from the nuclear encoded ITS regions including 5.8S 
rDNA of Dunaliella. Bootstrap values for 1000 replicates are given at the internal nodes. 



Indian isolates of D. viridis, while the presence of 2 stig- 
mata in S133 (only at lower salinity) has to be considered 
for re-examination and for final taxonomic identification. 

The remaining 2 strains S125 and S147 were placed in 
2 GROUPs (IV and V) as they were more clearly sepa- 
rated from other Dunaliella strains on the basis of 



genetic characters than morpho-physiological traits. 
Based on the taxonomic key [3] the strain SI 25 was 
identified as D. minuta (longer pyriform cells) but with 
clear separate starch granules in pyrenoid differing from 
D. minuta. This strain from Goa salt pan clustered with 
D. viridis (Figures 4 and 6), but with larger divergence 
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Figure 5 Phylogenetic tree (ML) of ITS region out grouped with D. tertiolecta; illustrated groups of 10 Indian Dunaliella strains. 

Bootstrap values for 1000 replicates are given at the internal nodes. 



values in both (ITS 4.98%, rbcL 3.84%) the phytogenies. 
Further due to the lack of molecular similarity with 
reported D. minuta (NCBI-BLAST analysis of ITS 2, 
results not given), the identity of the strain was kept in 
question and placed in GROUP IV. The identity of the 
Kutch strain S147 was a little confusing but interesting. 
It resembled D. tertiolecta in general morphology and 
18S rDNA size (-1820 bp), while grouped with D. viridis 
in ITS phytogeny (Figures 4 & 5) and with D. salina in 
rbcL phytogeny (Figure 6). It was isolated from a salt 
pan, having a dominant palmella stage and with a little 
higher /^-carotene content (6.7 pg.cell" 1 ) in stress, while 
D. tertiolecta was reported as a marine species without a 
palmella stage in its life cycle [3]. These observations 
resulted in the grouping of S147 separately as GROUP V 
and are showing a probability for a new species in the 
group. 

Diversity in Indian Dunaliella strains 

Buchheim et al. [9] have reported diverse community 
formation of Dunaliella in heteroclimatic hypersaline 
soils than in purely aquatic habitats. They hypothesized 
that external factors like temperature and salinity can 
enhance diversification and apparently got supporting 
results from the phylogenetic study of about 30 different 
isolates of Dunaliella (where 3 different morphotypes 
were characterized), based on 4 genes (18S, 26S, ITS & 
rbcL). Subsequently, Azua-Bustos et al. [40] reported a 
morphologically distinct, new Dunaliella species, D. ata- 
camensis, well adapted for sub-aerial life and with higher 
genetic divergence from its sister species. Our isolates 
are purely from aquatic habitats, but with high level of 
environmental fluctuations, especially in salt pans, and 



showed high divergence when compared to the reported 
Dunaliella species (Cladel of Figures 4 & 6) from NCBI. 
The geographic distance and isolation of the locations 
(Figure 1), from where the strains were obtained could 
be proposed as a reason for the divergence among the 
above Indian Dunaliella isolates. However 100% se- 
quence similarity and morphological resemblance 
observed between the two isolates SI 15 and SI 22 (iso- 
lated from Chennai - East coast and Goa - West coast 
respectively) need to be taken into account. 

Grouping pattern observed in the reported Dunaliella 
strains from NCBI in the cladistic studies ([27,29], 
present study) suggests a taxonomic revision of the 
strains especially when there are comments on confu- 
sion regarding the taxonomic status of many reported 
species. Consequently, Browitzka and Siva [3] have pro- 
posed an elaborate morphology/physiology based exam- 
ination of each strain in conjunction with molecular 
biology. However, among the 28 morphologically differ- 
entiated species [3,8], molecular aspects of only few im- 
portant ones have been extensively studied and reported, 
and a very large percentage still remains unexplored 
genetically. Hence, even after a detailed study based on 
morphology, physiology and molecular aspects, particu- 
larly to avoid misnaming, strain codes were assigned to 
our isolates which are more appropriate for comparative 
studies as well as for future communications. Morpho- 
logical and physiological study precisely groups 6 Indian 
strains into 2 sections - the carotenogenic Section 
Dunaliella (S089 and S135) and the non-carotenogenic 
Section Viridis (S115, S122, S133 and S125) [3,7]. The 
probability of the remaining 4 strains (S086, S118, S121 
and S147) to come under Section Dunaliella is much 
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Figure 6 Phylogenetic tree of the Maximum Likelihood (ML) analysis inferred from the rbcl plastid gene partial sequences. Bootstrap 
values for 1000 replicates are given at the internal nodes* show the position changes of S147 and S135 strains. 



higher as they are more carotenogenic (S147) and closer 
to D. salina (S086, SI 18, S121) in molecular analysis. A 
schematic representation of the diversity in Indian 
Dunaliella is given in Figure 7. 

The sequence diversity within the Indian Dunaliella 
strains was distinct when compared to the listed species 
of Dunaliella (ITS region & rbcL gene), and shows pos- 
sibility of the presence of multiple species in the group. 
Without the knowledge of sexual compatibility between 
the genotypes, it is not possible to determine whether 
this diversity is really representing a biological/evolu- 
tionary species or merely an intra-specific diversity (9). 
However, Coleman et al. [41] have demonstrated a con- 
currence between ITS sequences and mating ability in 



Dunaliella spp. [14]. Since, high level of sequence diver- 
gence observed among the Indian Dunaliella strains, 
could be correlated with sexual incompatibility, chances 
of more species/subspecies with respect to ITS phyl- 
ogeny seems to be a realistic possibility. Accordingly, the 
present study proposes the ITS region to be selected as 
a molecular marker in taxonomic delineation, which is 
smaller than 18S rRNA gene (with introns) and more 
diverged than rbcL gene. 

Conclusion 

The present study clearly shows high diversity within 
the Indian Dunaliella and reliability of 18S rDNA, ITS 
region and rbcL gene sequencing as a molecular tool in 
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Figure 7 Schematic representation of diversity of Indian Dunaliella. Grouping was done based on the morphology, 18S rDNA size variation 
(Figure 3) and ITS and rbcL gene phylogenies (Figures 4 & 5). 



species identification and genetic diversity studies. In 
a recent study, based on morphological parameters 
Jayapriyan et al [33] have denoted the presence of 5 spe- 
cies of Dunaliella (D. bioculata, D. tertiolecta, D. viridis, 
D. minuta and D. maritima) from India (east coast). 
However in the same study, 18S rDNA species specific 
fingerprinting using primers of Olmos et al. [31,38] have 
illustrated the same isolates as completely different spe- 
cies of Dunaliella (D. parva, D. bardawil and an uniden- 
tified Dunaliella sp.). Hence in the present study, for 
more clarity on the species lineages, along with morph- 
ology and 18S rDNA size, phylogenies based on a more 
diverse ITS region and a more conserved rbcL gene were 
also included, which otherwise are not available for In- 
dian Dunaliella. Consequently, presence of 5 or more 
species (or sub species), including 2 promising strains of 
D. salina (Section Dunaliella) and 2 D. viridis? (Section 
Viridis) strains, has got confirmed. The genetic 
characterization further helped in the separation of mor- 
phologically similar strains and in the clustering of In- 
dian strains of Dunaliella into 5 groups. Similarly, 
clustering of the reported species in a single clade (Clade 



1 with 100% similarity) in a both the phylogenies clearly 
emphasizes most careful recording of species names [3]. 
Hence, for resolving the issue prevailing in Dunaliella 
taxonomy and for elucidation of taxonomic species lin- 
eation of unknown Indian isolates, it is stressed to have 
further detailed molecular assessment coupled with add- 
itional examination of morphological (based on electron 
microscopy) and biological traits such as reproductive 
behavior (asexual- palmella, aplanospores etc.) and sex- 
ual compatibility. 

Materials and methods 

Sampling, isolation and culture conditions 

Water samples were collected from selected water bod- 
ies along the Indian coast during the months of February 
to July 2009 (Figure 1 & Table 1). Strains were isolated 
by serial dilution directly or after enrichment for a 
period of 1 or 2 weeks. Further purification was done 
with agar plating and picking the single colony to obtain 
unialgal cultures. Dunaliella salina CS265 was pur- 
chased from Collection of Living Microalgae, CSIRO, 
Australia and used as a reference strain. The strain S089 



Preetha et al. Aquatic Biosystems 2012, 8:27 
http://www.aquaticbiosystems.org/content/871/27 



Page 14 of 16 



(D. salinal, an old isolate from Chennai salt pan) was 
collected from Phytoplankton (aquaculture live feed) 
Culture Collection at CMFRI, Kochi. After purification, 
all the cultures were maintained in 75 ml modified 
Johnson (J/I) medium [10] in 100 ml flasks with 1.5 M 
(-200 ppt) NaCl, at temperature 25±1°C and at an ir- 
radiance of 40-50 umol photons m" 2 s" 1 supplied by cool 
white fluorescent lamps on a 12 h light: 12 h dark cycle. 
All cultures were sub-cultured once in a month basis. 

Microscopy and morphological study 

Live culture samples were examined using a Nikon 80i 
Research microscope (Nikon, Japan) with DIC (Differen- 
tial Interference Contrast) optics and images were cap- 
tured using Nikon DSFi le camera. Major taxonomic 
features observed include size, shape and colour of the 
cell, length of flagella, characteristics of stigma, pyrenoid 
and chloroplast and other cytoplasmic inclusions like 
refractile granules. Scalar measurements such as cell 
length and width, were taken from a minimum of 30 
cells from each strain randomly during mid growth 
phase immediately after fixing the cells with 1% Lugols 
iodine. The descriptive statistics such as minimum, max- 
imum, mean and standard deviation were estimated for 
the above scalar measurements. One way analysis of 
variance (ANOVA) was performed using SPSS (Version 
10.0) to identify whether there is any statistically signifi- 
cant difference among different Dunaliella strains for 
each character. 

Salinity tolerance study 

For salinity tolerance study, different Dunaliella strains 
were cultured in 5 salinity concentrations viz., 0.5, 1.5, 
2.5, 3.5 and 4.5 M NaCl in 150 ml (250 ml conical 
flasks) modified Johnson (J/I) medium. Other culture 
conditions like temperature and light were kept constant 
as given for normal culture maintenance. Cell character- 
istics like cell size and colour were examined at late 
growth phase under DIC microscope (Nikon, Japan). 
Cell count was taken on every third day using a Neu- 
bauer haemocytometer. Cell density was calculated and 
plotted against days of growth to obtain optimum salin- 
ity for each strain. 

/3-carotene analysis 

Beta carotene was estimated under normal (1.5M NaCl, 
irradiance of 40-50 umol photons m" 2 s" 1 ) and stressed 
(3.5M NaCl, irradiance of 100-150 umol photons m" 2 s" 1 ) 
growth conditions. Total pigment was extracted from 
4 ml culture at late growth phase (25 th day) in 4 ml ice 
cold 100% acetone. Liquid cultures were centrifuged 
(8000 rpm, 10 min.), the pellet washed with distilled 
water and re-suspended in ice cold acetone and left 
overnight at -20°C until the pellet became colourless. 



The extract was centrifuged at 5000 rpm for 5 min and 
absorbance was taken for the supernatant at 454 nm 
wavelength. Readings were compared with standard 
curve prepared with synthetic ^-carotene (Type 1, 
Sigma, USA) in 100% acetone as described by Hajezi 
et al. [42]. Cell density was calculated for the same day 
of extraction and /^-carotene was calculated per cell in 
picograms. 

DNA isolation 

DNA was isolated from 10 ml liquid culture at late 
growth phase following modified phenol-chloroform 
method of Wu et al. [43]. Cells were pelletized by centri- 
fugation at 6000 rpm for 5 minutes, washed in distilled 
water and re-suspended in 450 ul TEG (25 mM TrisHCl; 
10 mM EDTA; 50 mM glucose) buffer (pH 8) with Lyso- 
zyme (5 mg/ml) and vortexed with glass beads and then 
added 50 ul 10% SDS. The tubes were then incubated on 
ice for 10 min. and added 8 ul Proteinase K (20 mg/ml) 
and further incubated at 60°C for 60 min in a water 
bath. Once the cells were lysed completely, the DNA 
was purified following standard phenol/chloroform ex- 
traction and ethanol precipitation [44]. 

PCR amplification, sequencing and phylogeny 

A gene fragment of 18S rRNA was amplified using con- 
served primers MAI & MA2 [5,31]. Reactions were car- 
ried out in Veriti Thermal Cycler (Applied Biosystems, 
US) with a total volume of 25 ul containing PCR buffer at 
lx concentration with 1.5 mM MgCl 2 , 0.2 mM each 
dNTP, 1.5 U of Taq polymerase (Sigma, USA) 5 picomoles 
of each primer and 25 ng of genomic DNA. Thermal cyc- 
ling initiated with 3 min at 95°C and then 35 cycles of 30 
sec at 95°C, 30 sec at 52°C and 2 min at 72°C. Final exten- 
sion was for 10 min at 72°C. Amplified products were 
checked by electrophoresis in 1% agarose gel. The size 
(bp) of the amplified product was calculated by comparing 
it with standard molecular weight DNA marker (Step up 
100 bp DNA ladder, Merck, India) using the software 
Image Lab version 3 (Biorad, USA). 

Internal transcribed spacer (ITS) region (700 bp), in- 
cluding ITS1, 5.8 S rRNA and ITS2, was amplified using 
the primers ITS1 and ITS4 (17). Thermal cycling followed 
an initial denaturation for 3 min at 95°C and 35 cycles of 
30 sec at 95°C, 10 sec at 55°C and 45 sec at 72°C followed 
by a final extension at 72°C for 7 min. Partial (700 bp) re- 
gion of rbcL gene was amplified using the primers rbcL 
475-497 and rbcL 1181-1160 following Nozaki et al. [45] 
and Assuncao et al. [46]. The reaction mix composition 
for ITS and rbcL gene were the same as in the case of 18S 
rRNA gene amplification (given above). Amplified pro- 
ducts were tested on 1.5% agarose gel. 

All PCR products were purified using GenElute PCR 
Cleanup Kit (Sigma, USA) following manufacturers 
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instruction. Cycle sequencing was carried out using for- 
ward primers (MAI and rbcL 475-497) for 18S rRNA 
and rbcL genes respectively. Whereas the ITS region 
was sequenced using both forward (ITS1) and reverse 
(ITS4) primers. Sequences of DNA fragments were 
imported to BLAST [47] for similarity searches with 
available database at NCBI GenBank. The sequence was 
further aligned with the various available sequences 
(Table 4) of Dunaliella spp. and, Chlamydomonas rein- 
hardtii (ITS) and Paulschulzia pseudovolvox (rbcL) as 
out groups using the CLUSTAL-W algorithm [48] in 
Bioedit 7.0 (DNA Sequence Analysis Software package). 
To clearly illustrate grouping pattern in Indian Duna- 
liella isolates, D. tertiolecta was out grouped in ITS 
(Figure 5) phylogeny. Pair wise genetic distances among 
the different Dunaliella species and between the present 
isolates were calculated based on Kimura 2 parameter 
model for ITS region and Tamura 3 parameter for rbcL 
gene. The best nucleotide substitution model selection 
and phylogenetic analysis based on maximum likelihood 
(with 1000 boot strap replications) were carried out 
using MEGA 5 [49]. All the sequence information gener- 
ated in the present study were deposited in the NCBI 
database (Table 1). 
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